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INTRODUCTION 


The development of heavy transport and military aircraft with 
high performance characteristics has placed increased demand on 
landing gear and tire components. Specifically, aircraft tires 
have developed through the years to a point where the ratio of 
tire weight to aircraft weight has steadily decreased. In part 
this is due to improved tire cord materials, tire compounds and 
tire structural design, but in addition it has resulted from the 
use of aircraft tires closer and closer to their ultimate limit. 
This has resulted in a need for development of methods for pre- 
dicting the strength limit of aircraft tires on some rational 
basis other than indoor dynamometer testing. 

This problem of tire strength is not one such as is 
encountered in the strength of normal metallic structures. The 
reasons for this lie in the complexity of the material characteris- 
tics of the constituents in an aircraft tire. There are several 
factors at work here. First of all is the direct problem of 
loading of structure in the usual sense. In an aircraft tire 
almost all loads are carried directly by the textile cord struc- 
ture of the tire, commonly nylon as is presently practiced. There 
are some compressive loads carried by rubber elements in the 
tread but these are rather minor. In addition to the strength 
problem, aircraft tires generate substantial temperatures as they 
roll, both during takeoff and landing. The temperatures are not 
high in the usual sense of metallic materials, being of the order 
of hundreds of degrees Fahrenheit or less, but they are high 
enough to seriously degrade the strength of the load carrying 
textile components in the tire. Excessive temperatures can easily 
cause a tire to fail even if the stress levels are not overly high. 

In addition to the interaction of temperature and stress, 
even more complex phenomena are at work in the repeated daily 
use of aircraft tires. These phenomena are basically associated 
with the long term strength and adhesion degradation of an air- 
craft tire due to both stress and temperature, operating cycli- 
cally as the aircraft tire is used in a cyclic sense, as well 
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as the interaction of stress and temperature in such long term 
strength degradation. This is an area where not much hard infor- 
mation is available, and represents a region where the field 
usage of tires must depend on experience and judgment rather 
than upon science. Because of this strength degradation factor, 
the useful life of aircraft tires is limited, but these limits 
are not clearly defined and they vary from one aircraft to another 
and from one airline operator to another. 

At the present moment there are substantial advances being 
made in the tire industry in calculating stresses imposed on 
textile cord structures due to loading of a tire, such as in- 
flation and load carrying. However the area of tire temperature 
build up, which is the second major cause of tire failure, has 
been almost untouched and at the present moment there is no 
analytical technique suitable for predicting aircraft tire inter- 
nal temperatures on a rational basis. This means that the tiro 
designer or manufacturer is not able to compare the anticipated 
temperature build up in a new design with previous designs or 
with known temperature-strength characteristics of the material* 
which he wishes to use. This is a serious shortcoming in the 
design process, and is furthermore a serious shortcoming in the 
process of preliminary design of transport aircraft, where a 
knowledge of anticipated higher temperatures during normal take 
off and landing cycles would be an important asset in assessing 
the validity of proposed tire, aircraft weight and ground speed 
combinations . 

For these reasons, the present research effort is directed 
toward construction of a rational method for evaluating internal 
temperatures of an aircraft tire and for calculating these 
rapidly and easily so that the method may be used as a design l eel 
in both the tire and airframe industries. 


SUMMARY 


A method has been developed for approximating the internal 
temperature distribution in an aircraft tire due to free rolling 
under load in a straight line. The method is based on the concept 
of using an approximate stress analysis of a material point in the 
tire as it rolls through the contact patch and out into the region 
free of contact. From this approximate stress analysis the 
mechanical work done on this material point can be converted into 
the heat generated due to the loss characteristics of the materials 
at that point. This means that each material point of the tire 
will have its own cyclic heat generation each time it rolls through 
the contact patch and is stressed. 

Due to the high frequency of cyclic stressing compared with 
the relatively slow diffusion times of heat in polymeric materials, 
the temperature distribution in the tire is considered to be 
axisymmetric. 

With a rate of heat generation established at each material 
point on a cross section of the tire, !ieat diffusion throughout 
the carcass can be calculated by the usual numerical methods which 
will, of course, involve loss of heat in the tire through various 
mechanisms such as convection at the outer and inner surfaces, con- 
duction through the bead area into the flange of the wheel, and con- 
duction through the tread surface into the runway material. 

All of these mechanisms modify the normal diffusion process so that 
at any time, tire internal temperatures are dependent not only 
upon internal rates of heat generation but also on the appropriate 
loss characteristics at the tire surfaces. 

An extensive experimental program has been carried out Jointly 
between the University of Michigan and the Impact Dynamics Branch, 
Langley Research Center, NASA. This program consisted of tempera- 
ture measurements made at internal points in a number of 22 x 5.5 
RPR and 22 x 5.5 12PR aircraft tiros. Most of the experimental 
data was taken at the NASA Wallops Island Test Facility. This was 
done using the NASA tire test truck, and the tiros were operated on 
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a concrete runway surface. The University of Michigan data was 
taken on the 120 inch dynamometer drum at Wright Patterson AFB, 
Ohio, where higher tire loads could be obtained. 

All data was obtained by thermocouples buried in the tire 
carcass prior to retreading the tires. These thermocouples were 
connected to conventional recording apparatus, and gave continuous 
readings of tire temperature at selected internal points. 

Experimental data from these measurements is compared with 
calculation over a wide range of tire loads, inflation pressures, 
speeds and deflections. Generally speaking the analytical model 
proposed here gives good P£r • >oment of calculation with measured 
values when the tire is o;.„rating in the range of its normal rated 
load, pressure and deflection. When the tiro is seriously under- 
inflated the theory predicts lower temperature rise than is ob- 
served, and this aspect of it requires further work. 


GEOMETRY OF TIRE DEFORMATION 

The basic material property being used in the development of 
this computational method is the hysteretic loss characteristic 
of polymeric materials during cyclic stressing. .This is best 
illustrated in Figure 1, where a typical stress strain curve for 
a lossy material is illustrated. As material is stressed from 
point A to point B and return in Figure 1, the stress-strain 
curve on the loading cycle does not coincide with the stress- 
strain curve on the unloading cycle, and the net resulting area 
under the curve which remains as the material is returned to its 
original stress state, point A as illustrated In Figure 1, repre- 
sents a loss of mechanical work during the cycle from A to B and 



Fig. 1. Typical stress-strain curve of a material 
illustrating hysteretic loss 


return. This lost work is converted into heat and is either 
diffused through the element and eventually lost from its surfaces 
or must be accounted for in some way in a heat balance. From the 
analytical point of view, the work done in r stress cycle as 
experienced by a material point in passing through the contact 
\ patch is expressed most usefully by the use of a simple linear 

viscoelastic model such as the Kelvin-Voigt solid, and using such 
1 a model the energy loss per cycle may be written as in Equation 1. 

I AE = U v AV • tan6 (1) 


II 




5 


6 


where 

AE 

** Energy loss/cycle 


U 

Change in Strain Energy from 
A to B 


tan<5 

** Loss Tangent of Material 


AV 

* Volume Element 


This formulation presumes that the initial stress induced by 
tire inflation has no influence in a linear sense on the energy 
loss per unit stress cycle, and this may be shown directly by 
computation. However, the initia ? stress state may influence 
the material characteristics, which are not linear in this case. 
Methods for calculating the change in strain energy will be dis- 
cussed in the next section of this report. 

The fundamental method used for calculating the change in 
strain energy is based on two assumptions. The first is that the 
tire is made up of a series of material points, each of which 
undergoes a change in stress state corresponding to the points 
A and B in Figure 1, as it moves from point A to B shown as in 
Figure 2. 


A 



Fig. 2: Extreme points of stress excursion in a tire. 

As a material point moves from A to B in Figure- 2, we assume 
that the general wave form of the stress cycle depicted in Figure 
1 is unimportant, and we neglect higher harmonics of the cyclic 
stress in favor of a single fundamental cyclic stress change, 
which will be described in terms of an assumed geometry of the 
tire at point A and B. In other words, we make no allowance in 
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this simplified analysis for a Fourier decomposition of the stress 
state as the tire moves from A to B, but rather consider only the 
fundamental harmonic as contributing to the change in strain energy 
between the two points. This is felt to be a rational assumption 
in view of the fact that from physical considerations the funda- 
mental must be the primary component, and further that the inclusion 
of higher harmonics is not warranted within the framework of an 
approximate analysis such as to be carried out here. 

The second major assumption made in calculating the strain 
energy of the tire as it rolls through the contact patch is that 
at both points A and B of Figure 2, the cross section of the 
tire may be represented geometrically by its neutral axis, around 
which act both membrane and bending strains. This midline repre- 
sentation is assumed to be in the form of circular arcs, as shown 
in Figures 3A and 3B. The inflated shape of the tire is used to 

obtain best-fit values for the radii r and r , and for their 

J. c 

junction points 1 and 2 shown in Fig. 3a. 



(a) (b) 

Assumed geometry at * Assumed geometry at 

point A, Fig 2. point B, Fig. 2 


I, Fig. 3 

[’; 

!’ 

In constructing Figure 3 we illustrate the neutral axis of the 
structural part of the carcass by the lines labeled 0-1-2. The 
inflated tire geometry is presumed to be represented by Figure 
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3a, where a circular arc of radius r represents the tread re- 
gion, while the sidewall region is represented by another circu- 
lar arc of radius r^, tangent to the first arc at point 1. The 
deformed tire, under deflection A as shown in Figure 3b, is pre- 
sumed to operate under conditions of an indefinitely large radius 
in sections 0-1, and a reduced radius in the sidewall whose radius 
of curvature is given by the symbol r f . The bead, point 2, is 
fixed. The arcs 0-1 and 1-2 are allowed to shorten in length due 
to the* reduced tension in the sidewall. However, due to 
the geometry the shortening of 0-1 has little influence on the 
radius of arc 1-2, and for computational convenience the arc 0-1 
is considered inextensible. This allows the complete geometric 
description of the various radii of curvature, once the tire de- 
flection A and once the initial tire geometry are given under in- 
flated conditions. The computational details of the geometric 
representation are given in Appendix A. 

A material point is subject to both bending and membrane 
strain during deformation of the tire from point A to point B 
shown in Figure 2. The bending strains are dependent upon the 
radius of curvature of the midline of the tire carcass as illus- 
trated in Figure 3, as well as upon the distance measured per- 
pendicular from this carcass midline. Denoting this distance by 
the symbol y, the bending strain can be calculated in the usual 
linear elastic fashion based on a strength of materials approach, 
and based upon a knowledge of change in radius of curvature of the 
carcass midline between the two points in question in Figure 2. 

In addition, membrane strains exist in the tire cross section 
due to the high inflation pressure commonly used in aircraft tires. 
Some of these membrane strains are not directly obtainable from 
the geometric changes which have been postulated, since they depend 
heavily on contact patch boundary conditions. Other membrane 
strains must be obtained from a knowledge of the membrane stresses 
in the tire. These membrane stresses in turn depend upon the curva- 
ture characteristics of the tire, as well as upon certain geometric 
variables which may be calculated from the tire cross section geo- 
metries given in Figure 3. The details of these calculations are 
given in Appendix A, where a more complete geometric description 


9 


of the tire cross section is derived. In particular, the side 
wall area must be carefully defined since the point of vertical 
tangency of the tire is an important factor in computation of 
membrane stresses. This point of vertical tangency, shown as 
point 8 in Figure 3, shifts during the deformation process, and 
the position of the various material points relative to this 
point of vertical tangency must be tracked in the subsequent 
computations. 

Due to the fact that aircraft tire deflections are normally 
large, a further fundamental assumption is that the tire deforma- 
tion process causes changes in the tire cross section such as 
shown in Figure 3, such that deformations in the circumferential 
direction are negligible. This is equivalent to assuming that 
the tire is actually inflated in a barrel, so to speak, so that 
axisymmetry of tire displacement is maintained. A typical de- 
formation pattern illustrating this is shown in Figure 4. This 
imposes the necessity for tracking the radius of a material 
point from the axis of rotation of the wheel during the deforma- 
tion process, as shown in Figure 5. The variable used for such 
measurement is denoted py the symbol z. 
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Fig. 5: Tire Cross-Sections 

In Figure 5, point j represents a material point which changes 
axisymmetric radius during the deformation process from A to B. 
The details of the computation of such change are also given in 
Appendix A. 

Bead tensions also change during the deformation process 
and must be incorporated into a cyclic stres mechanism involving 
a loss characteristic of the bead itself. These bead tensions 
depend on angle of tangency of the side wall at the bead, denoted 
as point 2 in Figure 3. The details of computations of bead 
tension are given in Appendix B, where the bead geometry and 
assumptions governing it are discussed in more detail. 


TIRE DEFORMATION AND STRESS 

For the calculation of the strain energy change during the 
rotation of a material point through the contact patch we rely 
heavily on the assumed tire deformation given in Figures 2-5. 



Fig. 6: Tire Coordinate Notation 


Figure 6 illustrates the geometry involved and the notation being 
used here. R represents radial distance from the axis of ro- 
tation of the tire, <J> represents the meridional direction and 0 
represents circumferential direction. 

For purposes of computing changes in strain energy states 
during the cyclic loading process the tire is divided into five 
regions, denoted as I-V in Figure 7. 


I 


# 


Fig. 7: Tire cross-section showing regions of 

different characteristics and loading 
patterns . 

The rationale behind the division of the tire into these parts is 
that region II represents the tread region of the tire, which is 
considered to be rubber-like and in contact with the runway sur- 
face. It is usually characterized by material isotropy and by an 
elastic shear modulus which is nearly linear up to reasonably 
large strains. Further, it exhibits near incompressibility, and 
that fact is used in computing the strain state associated with 
contact against the runway surface. 

The region denoted by I is immediately beneath the tread 
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region of the tire and is firmly affixed to it, so that the 
strain state of the two elements in the membrane sense is 
considered to be equal. However, region I is assumed to be olane 
orthotropic in its characteristics, and may be thought of as a 
sheet of material lying in the <(>-0 plane as illustrated in Figure 


e 

w 
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Fig. 8: Laminate characteristics in Region I. 

In Figure 8 the diagonal lines represent textile reinforcing cords, 
which serve to give the primary stiffness to the carcass region. 
Because of their diagonal character, which is characteristic of a 
bias aircraft tire, the element shown in Figure 8 will exhibit 
different elastic constants in the <J> and 0 directions, as well 
as appropriate shear stiffnesses and Poisson's ratios which are 
well characterized by the mathematical model of a plane orthotropic 
material. Such materials have been widely studied in the tire 
industry and their characteristics can be estimated. Appendix C 
gives a detailed derivation of the appropriate elastic constants 
for a plane orthotropic material using the structure shown in 
Figure 8, adequate to within a range of diagonal reinforcing 
angles encountered in aircraft tire construction. Figure 8 can 
be characterized in terms of elastic constants dependent upon the 
cord reinforcing materials and upon the half-angle a illustrated 
in the figure. 

Regions III and IV of the tire, as illustrated in Figure 7, 


represent the side wall areas which are subjected also to changes 
of curvature as the tire deforms. These areas do not contact the 
runway surface, so that in region III the carcass strains are 
again determined and the elastic constants obtained from the model 
of Figure 8. These are utilized to determine the appropriate 
strain energy expressions. Region IV represents a rubber cover- 
ing used primarily as a scuff or abrasion resistant coating on 
the outside of the tire. It contributes to tire mechanical loss 
and to the subsequent heat build up, but is not structural in 
nature. We assume that it is bonded tightly to the carcass of 
Region III and undergoes the same membrane strains as those of 
Region III, but in addition undergoes bending strains due to its 
distance away from the neutral axis, which because of the rela- 
tively high textile modulus is usually located in Region III of 
the carcass. 

Finally, the beads of the tire undergo large tension force 
fluctuations, and these fluctuations must be accounted for separ- 
ately since they are caused by changing forces in the tire side- 
wall, which are reacted by the entire bead ring. This requires 
a separate analysis, not only for the force fluctuation in the 
bead ring but also because the bead rings are conventionally made 
of rubber coated steel wire which has significantly different 
elastic and loss characteristics from the remainder of the tire. 

For that reason this is denoted as Region V and is handled as a 
separate analysis. 

Because of limited computational ability, we approach the 
calculation of strain energy stored in the tire during the de- 
formation process by means of splitting each region of the tire 
into a number of cells, each cell being denoted by a pair of 
numbers in matrix notation fashion. This 'is illustrated in 
Figure 9, where cells are superimposed on the previous Figure 7. 

The elements shown in Figure 9 are purely for illustration purposes. 
Any number of elements may be chosen and the order of them may 
be adjusted to suit the needs of the user, but cor the instances 
to be illustrated numbering begins in the bead area and term* nates 
at the center line of the tire crown. Numbering of the vertical 
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location of the elements begins 



at the inner surface and terminates 
at the tread or outer surface of 
the tire, as the case may be. 

The inner liner of the tire, 
while pure rubber, is thin enough 
so that it may be lumped with the 
composite carcass material. 

Computation with the indivi- 
dual elements is facilitated by 
use of the matrix notation illus- 
trated in Figure 9. However, 
for a general theoretical treat- 
ment of the strain energy change 
in the tire in this report we 
will adopt the notation of Figure 
7 where the major portions of the 
tire are delineated. We 'now give 
an approximate analysis of each 
of the regions I-V of the tire. 


Fig. 9: Element numbering 

system 


Region I 

In this section we first used the assumed geometrical changes 
of Figure 3 to develop membrane strains in the <f> (meridional) and 
6(circumferential ) directions. We neglect strains perpendicular 
to the carcass midline. We include a shear effect to be discussed 
later . 

Consider first a side view of the tire as shown in Figure 10. 
Recalling the assumption of axisymmetric deformation, we completely 
restrain all circumferential motion so that the maximum membrane 
strain in the 0 direction occurs at point B and is given by 







A 



Fig-. 10: Side view of tire 

In the meridional, or direction the situation is more 
complex since we do not have access to convenient measures of 
deformation which lead directly to the strain. However the known 
curvature changes in the sidewall do lead to the membrane stress 
o., and in particular this stress can be evaluated at the edge of 
the contact patch such as shown in Figure 11 as points C and E. 

We now assume that roadway friction is not sufficient to 
substantially alter the value of across the relatively narrow 
contact patch width, and since the center of the contact patch 
is flat, then it holds that the value of o . at its edge must be 

<p 

constant across the entire contact patch. These membrane stresses 
are given later in the report as Eqn. (15). 



Fig. 11: Cross-section of tire. Points C and E are 
at edge of contact patch 
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In addition to these membrane strains, we assume that in 
region I the carcass becomes completely flat when in contact with 
the runway surface. This means that the curvature changes un- 
dergone by the carcass are exactly equal to the reciprocals of 
the principal radii of curvature of the tire, namely the values 
of R of Figure 10 and r of Figure 11. Using these values of 

C 

curvature change and the usual notation for bending strains, along 
with the modulii of elasticity for a plane orthotropic materials 
as illustrated in Figure 8, then the additional strain and the 

additional stress o. due to bending can be gotten and added to 

<P 

the e a and the o. associated with membrane effects. Note that 

when running on a test wheel, the drum curvature can be taken 

into account by defining the curvature change in the 0 direction 
1 2 

as (p+g). where D is drum diameter. 

The equations governing the deformed contact patch are now 
given by Hooke's law for a plane orthotropic material 


e 


0 


°e % 

E 0 " E <f> 
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(3b) 


where £g 
gives 


and o. are known. 
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Solving for the other two unknowns 


E e c 0 + e ^0#% 
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(4) 


and 


e ^ (1 ~ 11 <t»e e e 


(5) 
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These expressions now give the state of strain in region 1. 
Using the stress deviator and the strain state, the strain 
energy stored may now be written as 


U, 


M,B 


1 f 2 V* 

2]Ve + lj“ a “ w < 


+ °4, e © ( E^ 


r w 4>e ) 


( 6 ) 


While values of Eg may be obtained readily, the latter two terms 
of Eqn. (6) contain combinations of elastic constants which are 
difficult to evaluate. Net theory is of little help here since 
it predicts that both the terms 

and 

g 

<17 >W> < 7b > 

<f> 

are zero for inextensible cords. On the other hand, isotropic 
materials give a value of from 0.75 to nearly 1.0 for Eqn. (7a), 
but a value of zero for Eqn. (7b). 

Lacking better data on cord rubber composites, we must con- 
sider values of the constants given by equations (7) to be part 
of the numerical input data to any computational program. 

The other contribution to strain in Section I is associated 
with the fact that tensile forces in the region of the contact 
patch change fairly rapidly, thus generating shear stresses in 
the tire sidewalls associated with equilibrium considerations. 

This analysis is rather lengthy but leads to contributions which 
are in the present instance numerically small. In spite of this 
fact, they are included here since for tires of different geometric 
construction this effect may be substantial, and would be partic- 
ularly important in the thermal analysis of tires of radial con- 
struction. The details of this analysis are presented in Appendix 
D, but the resulting strain function is presented here in summary 
form; 


20 


U. 


r 2 

|[ 0 .857p o (l- ji)y s ] /(h 0 2 G„) 


( 8 ) 


where the symbols used are defined either in the notation given 
in the front of this report or in Appendix D, 

The net total strain energy change in Region I is now given 
by the sum of Equations 7 and 8 


U 


I 



+ 


(9) 


Region II 

In Region II the material is tread rubber, which is consider- 
ed to be isotropic and incompressible. For materials of this type 
it can be shown that the strain energy may be expressed as 

°II - G r[ E H +e y] • < 10 > 

where e., e Q> e. are the principal strains and G„ is the trend 
9 u y r 

rubber shear modulus. 

The membrane and bending contributions to the strains e^, 
and e « are exactly the same as those given in Equations 7 and 8 in 
the previous section, where now the variable y denotes the dis- 
tance of a particular computational cell or element above the 
neutral axis of bending of the section. In addition, the strain 
Ey. must be estimated from considerations of generalized Hooke's 
law as well as from contact stresses with the runway surface. 

To a first approximation, the runway contact pressures against 
the tread surface are close to the inflation pressure of the tire, 
particularly in the case of aircraft tires where inflation pres- 
sure effects are large compared to bending stiffnesses of the 
carcass. However it is known that the contact pressure distribu- 
tion in a tire on a flat surface is not uniform. Various experi- 
mental measurements have been made on actual tires under running 
conditions and these generally lend to the conclusion that the 
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shapes of the contact pressure distributions are approximately 
parabolic in the circumferential or 0 direction, while in the 
meridional direction, transverse to the direction of rotation, 
the distributions may range from approximately uniform, to para- 
bolic, to somewhat peaked in the shoulder areas. A numerical 
study of the influence of these various shapes on the strain 
energy stored in the elements in compression is carried out in 
Appendix E, where it is shown that to a first approximation 
a good numerical estimate of the influence of contact pressure 
on strain in the tread elements is given by: 


e 


y 



(ii) 


where E is the trend rubber modulus and p is the inflation 
r 

pressure. Combining Equation 10 with a generalized Hooke's law 
gives the total strain in a direction perpendicular to the carcass 
midline as Equation 11: 

i.5p 

E y - - -EJ-° - °-5<W <12> 

whore again and £g contain the membrane contributions from the 
carcass structure plus the bending contributions. Equations 11 
and 12 may now be combined into Equation 10 in order to form the 
total strain energy per unit volume associated with Region II of 
the tire. 

Region III 

In Region III of the tire we again consider the principal 
strains in the <p and 0 directions, but in this case there is no 
contact with the runway surface and strains in the direction 
perpendicular to the midplane are neglected. We include both 
membrane and bending effects in the calculation of the total 
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strain on each computational cell or element. The bending strain 
in the $ direction may be obtained from geometric considerations 
directly, since it arises from changes in radius of curvature 
which can be computed using the circular arc model discussed in 
Appendix A and illustrated in Figures 3 and 5. Bending strains 
in the 6 direction are considered negligible in this region. 

The membrane strain in the circumferential direction 6 is 
taken in this case directly from considerations of axisymmetric 
deformation of the tire. It is seen that as the tire passes 
from point A to point B of Figure 2 , the radial dimension of a 
material point j in Figure 5 changes from z^ to z^. This allows’ 
calculation of the membrane strain for an axisymmetric deformation 
pattern, as given by Eqn. 13. 


e e 


z r z i 

z i 


(13) 


The membrane stresses in the tire sidewall must be used to 
obtain the strains in the 4> direction. This is so because we 
have no ready geometric model for estimating these from the tire 
shape changes as illustrated in Figure 3. Referring to Figure 3, 
it may be shown from equilibrium considerations that the change 
in membrane stress in the direction between Figures 3a and 3fc 
may be written as 

. . p o 

% - Sir 

T \ c 


r f< 1+ z 


f b il 

f 1 z i 


(14) 


Since the strains in the 0 direction are known, Eqn. (13) and 
the stresses in the 4> direction, Eqn. (14), then Hooke's law 
must be used to obtain either the stress in the 0 direction or 
the strain in the $ direction. Again assuming orthotropic proper- 
ties for the sidewall structure, and solving for the strain in 
the $ direction gives once more, as in Eqn. (5). 
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C 4> 



(15) 


This leads to the sane expression for strain energy as was found 
previously for Region I, 

u m,b “ $j E e e e + E^ (1 “ M ^e p 64i )+o 4> G e ( E^ w 0<>“ y 4»e|| (16) 

In addition, the shear stresses caused by variation in membrane 
forces previously discussed as Eqn. 7 also appear in the sidewall 
in Region III. The expression for these is exactly the same as 
Eqn. 7, and will not be repeated here. However, using this term, 
the total strain energy change per unit volume in Region III may 
be obtained by summing Eqn. 16 and Eqn. 7 to give: 
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III 



+ u 


s 


(17) 


This implies that all of the shear forces generated by the vari- 
ation of membrane stresses are carried by the carcass region of 
the tire, Region III, and none by the sidewall rubber denoted as 
Region IV. 

Region IV 

In Region IV of the tire we again deal with a rubber covering 
but in this case one which does not contact the runway surface. 

The material is assumed to be isotropic and incompressible. The 
same considerations apply here as applied* in Region II except 
that the strain perpendicular to the carcass midline is neglected 
in this case. The modulus of this rubber covering is low enough 
so that it is extremely compliant as compared with the sidewall 
carcass structure. For this reason it is assumed that it will 
take on the same membrane strains as the sidewall carcass region 
III, and will only vary from it due to bending effects caused by its 
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further distance away from the neutral axis of bending. For this 
reason the very same expressions used for the calculation of the 
bending and membrane strains in Region III are also valid here. 
Specif ically, these are Eqns, 13 and 14, and using them the strain 
energy change per unit volume in this region may also be expressed 
in a fashion similar to that of Eqn. (10), again because of the 
fact that the $ and 6 directions are principal directions. This 
is given by Equation 18. 

B IV * G r[s 2+E e 2 ] < 1S > 

where is given by Eqn. 15 and by Eqn. 14, using the appropri- 
ate values of the variable y to denote the distance of the speci- 
fic computational cell or element from the neutral axis of bending, 
and further using the appropriate dimension z in Eqn. 13 for the 
same cell. 

Region V 

In this region we consider the bead elements themselves as 
components which lose energy during cyclic stressing of the tire. 

An analysis of this is given in Appendix B, and the result of that 
analysis may be used to introduce a total or gross strain energy 
change in each bead of the tire. Using the expressions derived 
in Appendix B, this becomes 

U v - 5< T Rf- T Hl> 2 /< A B 2 V 


(19) 


HEAT GENERATION RATES 


For purposes of obtaining heat generation rates from strain 
energy stored during elastic cycling of the tire we utilize 
the conventional Kelvin-Voigt mechanical model of an elastic 
solid which exhibits loss characteristics under cyclic stress. 

Ref raring to Figure 1, we see that it can represent experimental 
data observed on many materials, including the polymeric materials 
used in aircraft tires such as rubber, nylon in textile cord form 
and even the bead wire assembly, provided that the loss character- 
istics are measured properly. One widely accepted method of re- 
presenting such a phenomena as the Kelvin-Voigt model, is illus- 
trated in Figure 12 using the conventional spring and dashpot. no- 
tation. 



Fig. 12: Spring-dashpot model of a material 

exhibiting loss 

We impose on Figure 12 a cyclic sinusoidal strain in the form 

e s e sinwt (20) 

o 

The stress necessary to maintain such a sinusoidal strain is given 
by 


o * Ke + Cc “ Kc sinwt + Owe cosut (21) 

o o 

It is customary to convert this representation of Eqn. 21 into a 
form using the elastic constants E' and E", representing the char- 
acteristics of a material exhibiting loss, in the form of Eqn. 22. 
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o ■ E’e sinut + E"e coswt (22) 

o o 

It is clear from that representation that both E' and E" may be 
functions of frequency and usually are. In addition, particularly 
in polymeric materials, they show marked dependence on temperature 
and of course material composition. 

The energy loss during one cycle of stress is given by 

r 2 

Energy Loss ■ W ■ I ode » irE"e o (23) 

cycle 

We define the ratio of E" to E' as a material property called tanfi 
as given by Equation 24. 


E"/E ' = tanfi 


(24) 


Combining Equation 23 and 24 gives 

W ■ E'e o 2 tan6 ■ nUtanfi (25) 

where U is the elastic strain energy stored during the cyclic 
stressing process. This is the form in which energy loss is com- 
mon 1 v expressed, using the elastic constants E' and E” directly 
obtained from experimental data, or alternately using tan6 directly 
obtained from measurement. Several instruments are commercially 
available for the measurement of loss properties such as tan6 , and 
values of tan5 are commonly available for many polymeric materials. 

The materials used in tire construction exhibit properties 
which are close to being incompressible. For that reason it is 
assumed that the hydrostatic component of stress results in no 
elastic energy storage in a typical element of material, and hence 
exhibits no loss during the cyclic stressing process. Using this 
assumption, it may be shown that there are no effects of pre-stress 
on the energy loss per stress cycle, and for that reason Figure 1 
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of this report, showing the cyclic stress process between two 
general points, is drawn independent of the pre-stress level in 
the tire. In other words, we begin the strain measuring process 
from the inflated state of the tire and only consider that energy 
storage which occurs in a material element as it moves from the 
inflated but unloaded state to a position which lies directly 
in the center of the contact patch of the loaded region. The 
units of Eqn. 25 are the same as those of strain energy, namely 
work per unit volume. These must be converted to a rate of heat 
generation per unit volume in order to be used in the general 
heat diffusion equation discussed in the subsequent sections. 
Since aircraft velocities are commonly expressed in English 
units of feet per second, tan6 is dimensionless, and the tire 
diameter may conveniently be expressed in inches, and using the 
elastic strain energy in inch/pound units, then one may write: 


U/12 


ft. lbs, energy 
cu. in . volume 


/stress cycle 


1 3 

x yyg BTU/in /cycle 

3 

x 252 cal/in /cycle 

1 3 

x cal/cm /cycle 

(2.54; 


But = v Q where v q = aircraft velocity, fps 

r f = rolling radius of tire, inches 

: r o - 4 / 3 

a) = rad/sec 

12v o x 

- (r ' -A/3 ) x 2if cycles/sec. 


Loss/cycle = W = irUtanS 
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3 

Heat generation rate in cal/cm /sec * 

« * u x ts x m * (575? )3 11 3§ x cv^3) * ” tanS 

Uv tan6 Uv tanfi 

- ' 00988 (r°-A/ 3 y : °- 01 (J^S/3 

where r = outside radius of tire, inches, 
o 

A = tire deflection, inches 
v q = aircraft velocity, feet/second 
U = elastic energy, inch-pounds/cubic inch. 


( 26 ) 
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TEMPERATURE DISTRIBUTION MODEL 

Referring first to Figure 9, the tire must be divided into 
segments or cells in each of which the rate of heat generation 
is considered constant , and in which temperature will be calcu- 
lated on the basis of being averaged across the cell area. In 
assigning the geometric location of cells such as illustrated in 
Figure 9, we attempt to maintain material constancy across each 
cell, in the sense that the interface between the pure rubber 
region and the carcass region is used as a dividing line for a 
series of cells as well. This allows the properties of each cell 
to be essentially uniform, which is a considerable convenience 
in subsequent calculations. The stress state in each cell is 
calculated on the basis of the centroid of each cell area, and is 
considered to be averaged across the cell volume. The geometric 
location of the cell is taken to be the coordinates of the cell 
centroidal area such as shown in Figure 13. 



Fig. 13: Distribution of computational elements 

in the trend region 
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We have consistently utilized quadralateral elements in construc- 
ting the numerical analog to this solution. Thus each element 
either loses or gains heat from four faces, and this may be 
through conduction when the element in question is bounded by a 
similar element or it may be by convection from a surface face, 
or it may be by conduction to either the flange of the wheel or 
to the runway surface for those special elements which come in 
contact with them. 

For purposes of conduction to the wheel flange, a similar 
set of elements to those shown in Figure 13 must be laid out on 
the wheel flange using the appropriate thermal characteristics 
of the wheel material, so that heat flow can take place by direct 
conduction between the bead area of the tire and the flange. There 
will inevitably be a certain contact resistance between the bead 
seat and the flange material, but this factor is as yet unknown 
and is neglected in the subsequent computations. 

Contact with the runway surface is more difficult to quantify 
since it occurs only periodically. To a first approximation we 
treat that process in a similar fashion using the thermal con- 
ductivity of the runway material, but allowing only that fraction 
of the conduction process to occur as is proportional to the time 
spent by each surface element in contact with the runway. The 
remaining part of the cycle is occupied by those tread surface 
elements which contacted the runway surface now freely moving 
in the atmosphere, and during that fraction of the cycle the 
normal convection process from the surface occurs. 

Using the well known expression for diffusion of heat in a 
solid body under conditions of internal heat generation, such as 
given by Carslaw and Jaeger [1] , the equation for conduction of 
heat for a typical element is given by Eqn. 27. 

2 26 : 

KV e = o q (27) 


This may be written for finite difference purposes in the form 


At 


q + C ‘K-l,K,t (0 K-l,r $ K,H ) + <X K,K+l,)l (0 Ktl,r e K,^ ) 


+ Ct K,t-l,Jl (0 K,^l~ 0 K,£ 1 + “K.fc.f+l^K.i+l^K.P 


where the elements used In Eqn, (281 are illustrated in Fig. 14; 
and where 

- _ KA 

a m,n “ T 

with 


K = thermal conductivity 

A = area of contact between subscripted elements 
d = distance between centroids of subscripted elements 



Fig. 14: Notation for thermal elements 
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The basic method of computation is to use small time incre- 
ments and to calculate the temperature rise in eacn cell utilizing 
the temperatures of adjoining elements from the previous time 
increment computation. There are several factors which make this 
type of computation possible. First of all, it is well known 
in heat transfer calculations that if the quantity given by 
Equation 29 is large, then stability of such a time based forward 
computation is assured and convergence to correct results will 
be obtained. For that reason we have chosen to keep cell elements 
moderately large and time increments small since this reduces re- 
quirements for computer memory, and also insures stability of the 
computation . 




M = 

Ax 



aAt 

where 

a 

- 4/« 



Ax 

= element size 


At 

- time 

increment 


Secondly, we are dealing here with a type of computation 
where the dominant term is the rate of heat generation q, and for 
purposes for aircraft taxi and take-off the heat build up occurs 
over a relatively short period of time, so that the heat genera- 
tion term almost completely controls the temperature rise in the 
tire. Heat diffusion, at least for the take-off case, is a mini- 
mal effect. This means that heat conduction using adjoining cell 
temperatures from the previous time interval incur only very small 
errors in the final temperature distribution in the tire. 

The heat transfer coefficients from tire to air are not 
particularly well known in this case, but again it is fortunate 
that taxi take-off cycles are relatively short, so that again the 
rate of internal heat generation is the dominant term. However, 
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in studying taxi and combined taxi take-off operations, more 
emphasis is placed on the accuracy of the heat transfer and heat 
conduction coefficients associated with this problem. In that 
case those coefficients become more important. The numerical 
proportioning of the heat transfer coefficients is given in detail 
in Appendix H. 

The method of temperature computation using Eqn . 27 now 
consists of starting from whatever initial temperature state one 
wishes in each of the cells, allowing a short time interval to 
take place and computing the heat generated in each cell using the 
concepts developed in the previous two sections. In carrying out 
such computations a knowledge of the material loss characteristic 
tanfi must be known, and these are given in the last section of 
this report. 

At the end of the first time interval the temperature field 
in the entire tire is calculated and is used as a basis for the 
heat transfer characteristics during the second time interval, 
which follows immediately. The computation is sequential, thus 
requiring minimum memory storage in a computer, and may conveniently 
be carried out on relatively unsophisticated equipment such as 
microcomputers or minicomputers. 

The contained air temperature in the tire is often measured 
and used as an indication of overall tire heating. It's calcu- 
lation requires a separate analysis given in Appendix H. From 
this analysis, the contained air temperature is shown to lag 
behind the inner surface temperature of the tire. In the initial 
stages of tire heating, its value is given by 




Ay = 0At(0-y) 

where 


• 


y 

= contained air temperature 


0 

*» weighted average of tire inner surface 
temperatures 


At 

= time increment 


0 

= ° constant involving properties of the 

air and heat transfer coefficients. 
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The heat transfer coefficients chosen for the inner surface 
are presently based on the concepts that inner peristaltic pump- 
ing velocities are small, just sufficient to insure good mixing. 
This is an area where more measurements are needed. 


EXPERIMENTAL PROGRAM 


The problem in hand is an extremely complex one and could not 
be expected to be verified without considerable experimental data. 

A number of assumptions are used in the computation of the rate 
of internal heat generation in the tire, and many of these assump- 
tions can only be verified by direct measurement of temperature 
fields in actual tires. For that reason a comprehensive experi- 
mental program was carried out by the Impact Dynamics Branch, 
Langley Research Center, National Aeronautics and Space Admin- 
istration, Hampton, Virginia, Most of the work was done at their 
Wallops Island Test Station, using a specially designed truck 
with a yoke and loading apparatus designed for the testing of 
tires . 

Thermocouples were installed at the University of Michigan 
in approximately twenty 22 x 5,5 8-ply and 12-ply rated aircraft 
tires furnished by NASA Langley, Most of these thermocouple 
installations were achieved by implanting appropriate thermo- 
couples in holes drilled in the tire carcasses after having been 
buffed but prior to retreading. Subsequently the tires were re- 
treaded by a commercial retreader, and were wired for internal 
temperature measurements during both straight line rolling, cor- 
nering and braking. For purposes of this report only the straight 
line rolling data will be used. Due to load limitations of the 
test equipment the full rated loads of the tires could not be 
applied, and since it was desired to achieve normal tire 
deflections, inflation pressures were also kept below normal values. 

The tire temperature data obtained from typical test runs is 
presented in the next section of this report where it is compared 
with computation. In general the use of* appropriate thermocouples 
carefully implanted in the tire was sufficient to give a good 
picture of the internal temperature distribution. A thorough 
summary of this information has been presented in a separate 
NASA document. 

Due to limited load and pressure combinations which could be 
used on the NASA equipment, the University of Michigan conducted 
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a short series of taxi-takeoff roll tests on ins t rumen ter* tires 
of the same size at Wright “Patterson Air Force Base, Dayton, 

Ohio, using their 120 inch aircraft tire dynamometer. With this 
equipment full rated loads and speeds could be obtained, This 
data is also presented in the subsequent section where it is com- 
pared with experiments, and this represents a different set of 
temperature rise values since it involves loads, pressures and 
speeds which are consistent with real operating conditions. 

Referring to earlier sections of the report, it is clear 
that the rate of heat generation depends on the strain energy 
stored in the tire, which in turn depends both on bending and 
membrane effects. While deflection itself controls the bending 
contribution to this problem, the membrane contribution is clearly 
dependent on the inflation and thus both effects are needed to 
give a realistic description of tire temperature rise during ser- 
vice conditions. 

While a great deal of temperature data was obtained by these 
two series of test programs, it is so voluminous that only selected 
parts of it have been used for purposes of comparison with calcu- 
lation in this report. In order to provide a complete record of 
the test data available, however. Table 1 shows the load and 
speed combinations available in straight line rolling from the 
NASA Wallops Island tests, while Table 2 shows the load, pressure 
and speed combinations available from the University of Michigan 
tests carried out at Wright- Patterson Air Force Base. Some of 
this data will now be extracted in detail and used in a subsequent 
section of this report for comparison with calculation. 



Table 1 

Operating Conditions NASA Wallops Island Tests 



Load 

lb 

Pressure 

psi 

Speed 

MPH 


2700 I 

83 

20 


3270 

83 

20 

NASA 

4000 

83 

20 

8PR 

4000 

115 

20 


4000 

150 

20 


2560 

55 

20 


3150 

55 

20 

NASA 

4000 

55 

20 

12PR 

4000 

110 

20 


4000 

78 

20 


4000 

78 

35 


4000 

78 

50 


Table 2 

Operating Conditions - University of Michigan Tests 



Load 

lb. 

Pressure 

psi 

Speed 

MPII 



4350 

125 

20 


UM 

3263 

87 

20 


8PR 

4350 

156 

20 



4350 

94 

20 



UM 

12PR 


7100 

5325 

7100 


209 

132 

261 


20 

20 

20 










CALCULATION VS, EXPERIMENT 


We now desire to compare calculation of temperature distribu- 
tion in the tire with data measured from both NASA Wallops Island 
tests and the University of Michigan tests. In order to carry 
out the computations the material characteristics must be assigned. 
For this purpose we use data estimated from the literature for 
the material characteristics of the tire as given in Table 3. 

Using a time increment, as defined in Eqn. 28, of t ■ 10 se- 
conds, and using the material characteristics and tire geometries 
for both the 8 PR and 12PR 22 x 5,5 tires, as determined in Appendix 
H, computation was carried out oyer short time intervals for the 
temperature rise in these tires under conditions experienced by 
both the Wallops Island tests and the tests conducted by the 
University of Michigan, Comparison between calculation and ex- 
periment is shown in Figures 15 through 33, where the experimental 
data is plotted as a function of time. The computations are 
carried out using Eqn. 28 and the material constants given in 
Table 3, along with the specific g^cmetry for these tires, as 
defined in Appendix H, 

It should be emphasized that while the outside dimensions of 
these two ply-rating tires are the same, their internal character- 
istics are quite different and the modeling is completely differ- 
ent. The elements representing the heat transfer cells are not 
the same, nor is the nature and extent of the carcass area. The 
8PR version of the tire has only one bead, while the 12PR version 
has two beads. For these reasons the two sets of comparisons 
should be considered as being taken on two separate tires. 

As can be seen from these figures, the calculations generally 
agree well with the experimental data. In particular they predict 
the ranking of hottest areas in the tire at a particular location, 
i.e. from outside to inside and vice-versa. They support the re- 
versal of inside-outside temperature gradient between sidewall 
and tread, a phenomenon observed experimentally. In addition 
they generally support the experimental data observed over a rather 
wide range of inflation pressures and tire deflections. 
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Element 

E 

G 

m 

a 

K(3) 

(4) 

P(3) 

Cp(3) 

Tread(II) 


335 

0.2 

- 

5xl0" 4 

1.0 

m 

Carcass 

(1) 


0.1 

(.2) 

5x1 0“ 4 

1.0 

■H 

(I, III) 







mm 

Sidewall 


335 

0.15 



1.0 

iWi 

Rubber CIV) 







■HI 

Bead (V) 

15xl0 6 


0.03 



" ■ -*-r 

■1 


Material Characteristics 


Table 3 


♦ 

See Fig. 7 


a) 

See Appendix C. Use G * 335. 


(2) 

Calculated by cosine law using o * 35 

0 crown angle 

(3) 

See Appendix H 


(4) 

cal-cm/°c-cm^-sec 


(5) 

Values of tan6 # are given for 25*C. 

For higher temperature 


tan6 is calculated by an approximate 

expression 


tan5 = tan6» e 
where AO is the temperature rise in 


- . 01 A0 

*C above ambient. 
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From a study of the comparisions between these computations 
and the experimental data it is concluded that the major trends 
in the temperature rise problem can be calculated adequately for 
relatively short term taxi*takeoff conditions. More emphasis 
must be placed on obtaining good thermal properties of the mat* 
erials and better heat transfer coefficients before we can hope 
to achieve success in calculating long term or near equilibrium 
temperature conditions in the tire. 



* 


Fig. 15: Measured and Calibrated Tire Temperature Profiles. 




asured and Calibrated Tire Temperature Profil 







Fig. 17: Measured and Calibrated Tire Temperature Profiles. 




Measured and Calibrated Tire Temperature 









Measured and Calibrated Tire Temperature 


























TIRE: 22x5.5-8PR ° INSIDE SURFACE: CONCRETE RUNWAY 

SPEED’ 20 MPH A MIDLINE Fz : 4000/4350 = 92% 

FREE ROLLING □ OUTSIDE 8 Z : 1.050/1.344=78% 



Measured and Calibrated Tire Temperature Profiles. 







TIRE : 22 x 5.5 12 PR O INSIDE SURFACE : 120 IN STEEL DRUM 

SPEE0:20MPH A MIDLINE F z :7100/7100 100% 

FREE ROLLING □ OUTSIDE 8 Z * 1.304/1.312 99% 

UM TEST CALCULATED Po * 209/190 110% 
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Measured and Calibrated Tire Temperature Profiles. 

















!2x5.5-!2PR O INSIDE SURFACE: CONCRETE RUNWAY 

20MPH A MIDLINE F z : 2560/7100 = 36% 

ROLLING O OUTSIDE 8z : I.C25/I.3I2 = 78% 
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Measured and Calibrated Tire Temperature Profiles 






Fig. 28: Measured and Calibrated Tire Temperature Profiles. 















22x5.5 -!2PR O INSIDE SURFACE: CONCRETE RUNWAY 

20MPH A MIDLINE F z : 4000/7100 = 56% 

ROLLING □ OUTSIDE 8 Z : 1.435/ 1.312 = 109% 

TEST — CALCULATED P ft : 55 / 190 = 29 % 



Measured and Calibrated Tire Temperature Profiles 









Fig. 30: Measured and Calibrated Tire Temperature Profiles. 


TIRE: 22x5.5-12 PR O INSIDE SURFACE: CONCRE r E RUNWAY 

SPEED: 20 MPH A MIDLINE F z : 4000/7100 = 56% 



Measured and Calibrated Tire Temperature Profiles. 













SURFACE * CONCRETE RUNWAY 



mm 


Measured and Calibrated Tire Temperature Profiles. 















Measured and Calibrated Tire Temperature Profiles 
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APPENDIX A 


Idealized Cross-Section Curvatures 

It is assumed that the geometry of the tire carcass can be 

described by two intersecting circular arcs. For the undeformed 

tire these arcs describe the carcass centerline from the crown 

of the tire to the shoulder, R , and from the shoulder to the 

c 

bead, R^, as shown in Figure Al. 

It is assumed that the x-y coordinates of the shoulder 

(x ,y ) and bead (x n ,y n ), the y-coordinate of the carcass at 

the crown, y , and the length of the carcass centerline from 

the shoulder to the bead, S T , are known. All of these geometric 

properties can be measured from appropriate tire cross-section 

profiles. The location of the shoulder is somewhat arbitrary. 

However, for this analysis x was taken to be 0.75 of the section 

c 

half-width w, and y to be the intersection of x=x and the 
carcass centerline. 

From geometry: 


*c * W 

2(y -y ) 


-1 x c 

e c = sin ( r^ 


s„ - R *8„ 

c c c 


(Al) 

(A2) 

(A3) 


To find the radius of curvature of the undeflected side- 
wall one must first determine H: 

H = Ax c -x B ) 2 +(y c -y B ) 2 (A4) 


Noting that = Sj/0 and that H/2 = sin (-g) , the angle 0 can 
be found from: 
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a 

6 - sin(|) = 0 ( A5) 

Solving (A5) for 0 then allows the calculation of R^. 

To calculate the radius of the deflected sidewall it 
is first assumed that the deflection is greater than y -y 

IT) C 

and that the carcass changes from a circular arc of R to a 

flat surface of half-length S . It is thus assumed that the 

x-coordinate of the shoulder in the deflected position is 

S and the y-coordinate is (y -A). With the new co-ordinates 
c m 

of x„ and y a new H can be calculated from A4. However, 
c c 

jefore a final radius of curvature of the sidewall can be 
determined, one must account for the approximate change in the 
arc length S^. This change takes place because there is a 
significant strain in the meridianal direction for the operating 
conditions usually encountered in aircraft tire service. 

Referring to Figure A2, the final length of the carcass 
centerline is: 

s t - s l< i+ v 

where is the meridianal strain. This strain is approximately: 


c <p 


Thus: 

P n 

S f = S L (1 + h-f-CW) ( A6 ) 

c <{> 

where P Q is the internal pressure and E^h c is the average stiffness 
of the sidewall in the meridianal direction. Also from the figure: 


0 f 

R f sin(^) 



and 



Substituting from (A6) and (A8) into (A7): 


(A7) 


(AS) 



Figure A2 - Geometry of original and final idealized curvatures 
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Rj sin 


s lU + 
[ 


idH R r R i» 

$ c 


] 


2Rj 



(A9) 


The solution of this transcendental equation for R^ gives the 
radius of curvature of the deflected sidewall. 


APPENDIX B 


Bead Tension Analysis 

Consider the sidewall of an aircraft tire as shown in Fig. 
Bl, first in the undeformed and then in the deformed state. 


<L <t 



Figure Bl . Sidewall geometry used in bead tension analysis. 
The membrane forces in the 4> direction are given in Eq. 

( 15 ). 


The angle p 

1 

Fig. Bl 


N 


♦ 


i 


N 


♦ 


f 


is given 


P 0 b. 

= s 2 ^ 1 + 4 >] 

P„ b * 

- ^V 1 + ^ ] 

from initial tire geometry. 


From 
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2r^ cos - 2r f cos B f 


A ' 


(Bl) 


or 


cos B, 


2r^ cos Bj - A 

2FZ 


(B2) 


Also 

A ' » A - r (1-cos X) (B3) 

c 

Using geometry previously developed in Appendix A, a' can 

be found from Eq. (B3), and B f from Eq. (B2). 

Bead tension may be split into radial aci axial components 

as in Fig. B ?. The axial component T is carried by the wheel 

z 

flange. The radial component causes bead tension in the amount 

T r *= (N^ sinB) R b (B4 ) 

using a uniformly loaded ring 

The difference in bead tensions 
between points A and B of Fig. 

2 is given by 

Fig. B2 

T R f - '''Rj ■ <% Si " «f - \ si " 8 i> R B (B5) 

Denote the gross area of the bead wires as A D , and the 

O 

modulus of the bead bundle based on this gross area as K_. Then 

B 

the strain energy stored during one half cycle is 

u v ” 2 (t r " t r. , 2 /(a b k h ) 

f 1 



(B6) 


APPENDIX C 


Aircraft tire construction is entirely bias ply and the 
ply angles span a limited range, say between 30° and 60° for 
half angles a as illustrated in Fig. 8. Previous studies by 
the authors [2,3] and by Walter [4] have shown that within this 
range of angles the elastic constants of the composite bias 
ply structure are well represented by a simplified analytical 
Tiodel using cords which are inextensible but which are encased 
in an elastic matrix. 

Consider the bias cord structure of Fig. Cl. 

I X Assuming inextensible cord, 

x ■ l cos a 
dx c -i sin ada 

= ^ = -tan ada 

♦ 

= i sina; 

■ l cos ada 
“ ^ * cot ada 

Cl 

l i 1 of idealized bias net structure. 
♦ 2 

Cy tan ac„ or 
2 

tan e^ 




c 6 4 2 
— = tan a = 


0 4 > 


cot a = 


te 


(Cl) 

(C2) 


OR 


Calculate the modulus of elasticity, assuming energy storage 
by rubber only. 



V t0 " 



Fig. C2 . 

■<j> Mohr's circle 

for idealised net 
structure. 


The radius of Mohr's circle is e. ( 1+tan a ) = (y/2). The center 

n T 

is located at ( 1-tan q )e^ = e 

The element is now subjected to both shear y and extension e. 
Express the strain energy as 

„ 1 ^ 1 ^1 
V = o e a +t ? e a + oY t 

2 x x 2 y y 2 'xy xy 


But t = Gy and also due to isotropy 
- xy xy 


cr v =cr 

i 


o 0.5a 
y E = — E" 



a = 2Ee 


Cy = a But G = 


E 

2(l+n) 


• J 


Th< strain energy is now 
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V = | • 6Ge 2 + | 6Ge 2 + ^ Gy 2 


sa 


^[l2G»( 1 ~ t | 1 > --) 2 + G( l+tan 2 a) 2 ]e 2 




4G[l-tan 2 a+tan 4 a 


3 ® 


2 

♦ 


E. * 4G[l-tan 2 a + tan 4 a] (C3) 

Also tan(90-a)=cot ct; Eg = 4G[l-cot 2 a + cot 4a] (C4) 

This agrees well with more exact expressions of Refs. C2-4D. 

The angle distribution in a bias ply aircraft tire is given 
by 

cos a - — cos a 
r o 



where 

= crown angle 

o 

r = crown radius 
o 



APPENDIX D 


Sidewall Membrane Shear 


The membrane forces and N Q vary as the tire rolls through 
the contact patch. The variation in Ng causes a shear stress 
to be required for equilibrium purposes. This shear stress gives 
rise to an energy storage, hysteretic loss and heat generation 
process which is similar to those found by other mechanisms. 

Consider the tire cut at its centerline and spread out 

flat as shown in Fig. D1 . Equilibrium of forces requires the 

existence of a membrane 

shear force N , as shown, 
xy 

, This is caused by the 

/N-ds 

t o variation in N_ with circum- 

9 

ferential position. 


/(N^+ANpJds 



Fig. Dl, Idealized view of tire section. 


Next consider the relation between N and N . For a net 

e <t> 

structure in which loads are carried by the cords alone, Fig. D2 

illustrates the division of cord load between N, and N„ . If T 

4 > 0 

represents a typical cord load, it may be shown that 



Fig. D2 . I d ea 1 i z od net 
structure . 
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Since is calculated as part of the contribution to strain 

energy in the sidewall regions I and III, then it is available 

for computation of N«. At any location the variation of N a 

o o 

follows the variation of N^, which is approximately of the form 
shown in Fig. D3. 



Letting the variation of with 
respect to circumferential position 
x be approximated by 


N e 


a + b(l-e~ c 



(D2) 


Fig. D3. Variation of N 0 

with circumferential position. 


* 2 

where a = (N^) B cot a 


b =■ [<V*-(V B ]cot 2 a 


2 1 

c 2 — ~ where w = section width, 
w^ 

Since the maximum value of N X y depends on the maximum 
gradient of Nq , then it can be determined by differentiation 
of Eq . D2 . 


dN 

dx 


« 2 2 

0 0 . 2 -c x 

- = 2bc xe 


d n ? 2 

— ^ = 2bc 2 tl-2c 2 x 2 ]e c x 

dx 


The maximum of dN occurs at x 

e 

dx 


1 

c/2 


♦See Fig. 2 for location of points A and B. 
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and results in 


dN fi 

< a3T ) max 



0.858 - 
w 


0,858 

w 


C(N 0 ) a -(N 0 ) b ] cot 2 a 


From plane equilibrium, 

dNg 

dx 



so that 

dN xy * C (N ^A- (N »>b] cot2 whe ‘’ e 

is measured along the circumference of the tire meridian. 

varies with this position, so that the net value of N at 

xy 

computational cell is given by 


N 


xy 


0.858 

w 




cot a-dy 


From this the strain energy can be formed 


- 1 »j 2 /, 2 — 

2 N xy/ 1 ^xy 


(D3) 

(D4) 

y 

Cot a 
any 

(D5) 


U 


(D6) 


APPENDIX E 


Contact Pressure Analysis 

Independent of other effects, the tire cycles between 
zero tread compressive stress and a maximum compressive 
stress. 

Measurements indicate a maximum pressure/ unit length at 
the b of contact, approximating 1.25 x average inflation pressure 
x width. 

Examine three possible pressure distributions, all giving 
rise to the same load/unit length. See Fig. E-l . 


Stored energy ~ 


w/2 


P 2 dy 


-w/2 


assuming E = constant across the 
tread . 


(a) E 


f w 

a - j 0 (1 — -o 


25p ) dy = 1.5625P w for constant pressure 


(El) 


(b) Parabola 1; p - p m [l-(|^) 2 ] 


But 


f w /2 2 2 r 

I p^dy = prn (y - 

8y 3 + 16y 5 
3w 2 5w^ 

w/2 

-w/2 


-w/2 

f w/2 

i pdy = 1.25 p w = 

rw/2 

j PJ 1 - ' 

i 2 

-w/2 

-w/2 

vv 


2p w 
m 


74 


75 


P m - | x 1.25 

m £ o 

E b ■ 4 * T5" - ! * § Po * fr " 1875 Po w ‘ 


(c) Parabola 2: p 

fW/2 

I pdy « 1.25 p Q W 

-w/2 


m 



m 


\ Sr dy 


m 



w/2 

-w/2 


2 , 

P dy 


= P, 


m 


16y 


w 


dy = 


2 w 
p m 5 


2.8125 p w 
o 


We adopt a conservative value of — - — — = 1 . 5 
case. Assume uniform contact pressure of 1.25 x p Q . 
multiply by 1.5 factor for shape distribution. 


(E3) 

1.2 x E o 
a 

(E4) 

= 1.80 E 

a 

x uniform 


Then 



APPENDIX F 


Meridianal Membrane Strains 


Consider the meridianal membrane stresses in a torus of 
circular cross section as shown in Fig. FI, subject to an 
internal pressure p. Cut out a section bounded by a vertical 
tangent point, distance b from the axis L, and a general point 
R. This is shown in Fig. F2. 



Horizontal equilibrium of Fig. F2, considering the entire 
asixymmetric ring, lends to the expression. 


N^* 2tt R sin $ = p.7t(R 2 -b 2 ) 


(FI) 


from which 


N 


p(R+b)(R-b) 
2R s in 4> 


p(R+b)(rsin<J)) 
2 R s i n <ti 


♦ 1 > 


(F2) 
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During deformation of the tire into the center of the contact 
patch, point B of Fig. 2, the tire deflects but the sidewall 
remains a circular arc, although of reduced radius. The vertical 
tangent point B also shifts inward, so that both r and b in 
Eq. F2 take on new values. The final value of N still carries 
the same form, namely 




(F3) 


where r f , b f and R f define the same material point previously 
denoted by r, b and R. The change in now becomes 

iN. - N. -N. (F4) 


This is the stress change upon which the strain energy is based. 


APPENDIX 0 


Contained Air Temperature 

Consider the tire cross-section as shown in Fig. G1 . Let 
0 be the area weighted average of the temperatures on the inner 
surface of the tire over the region 0-1-2. Let y be the contained 

air temperature. Then the equation 
for temperature rise of the contained 
air is 

V p ‘ v Ht ■ hA < 0 -y) (on 


Fig. Gl. 
where 

h = average heat transfer coefficient 
A = surface area of convection 0-1-2 
0 ^ = specific heat of air 

p = air density 
V = volume of contained air 

Let hA/CppV = 0. Then Eq . Gl becomes 

+ By = BO (G2 ) 

In the initial stages of tire heating the tire inner surface 
temperatures rise almost linearly, so that 0 = ct. In this case 
Eq. (G2) has a solution 


1 



y 


0 - |(l-e" Bt ) 


(G3) 
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This solution lags behind the tire inner surface temperature. 

For longer times when 9 is determined numerically, the 
contained air temperature y must be calculated by the methods 
previously described using a finite difference form of Eq. (Gl). 


APPENDIX H 


Geometry of 22 x 5.5 Tire 

Figures HI and H2 show the 22 x 5.5 OPR and the 22 x 5.5 12PR tires 
in cross-section used in the thermal analysis discussed below. 

For this analysis it is convenient to divide the half cross-section 
into eleven sections in which each section is in turn divided into 
three elements. As can be seen from the figure, each element is 
identified by two subscripts, the first one indicatirj the merid- 
ianal location from bead to crown and the second indicating the 
location through the thickness from inside to outside. Many 
geometric values are required from this sectional profile in the 
analysis which follows. The definitions of such values are listed 
below and their numerical values for both tires are given in 
Tables HI and H2. 


V(I,J) 
LH( I , J ) 

LV ( I , J ) 

DH ( I , J ) 

DV(I,J) 

AA(I) 

DA(I) 

AC( I ) 

DC( I ) 

R(I) 

YS( I , J) 
Y ( I , ) 


mm 

H2( I ) 


Volume of element 

Area across which heat is transferred in thickness 
direction within the carcass 

Area across which heat is transferred in the meridianal 
direction within the carcass. 

Distance between nodes through thickness 

Distance between nodes meridianal direction 

Area across which heat is transferred to outside air 

Distance fron, outside element nodes to outside surface 

Area across which heat is transferred to contained air 

Distance from inside element nodes to inside surface 

Radial location of each section from wheel centerline 

Distance along meridian from .crown to node of each element 

Distance from inside surface of section to node of each 

element in the section 

Thickness of carcass in each section 

Thickness of rubber in each section 



Figure HI. Element sectioning for Thermal Analysis 
Solid dots represent approximate thermo 
couple locations. 



Figure H2 - Element Sectioning for Thermal Analysis 

Solid Dots Represent Approximate Location 
of Thermocouples. 22 x‘ 5.5 12PR 













































APPENDIX I 





Thermal Analysis 


In the thermal analysis discussed below account is taken of 
conduction to neighboring elements, conduction to the road 
surface, convection to the outside air and to the inside air, 
and heat generation within each element. In writing the heat 
balance equation for each element it is necessary to provide 
effective heat transfer and heat generation coefficients for each 
element. The necessary coefficients are: 


K p * a • At 

AV, AH = — 

C R' t V v ' d 


QB 


q • v • At 
C R ,p R* V 


HR 


1 


1 

V a 

^PR^d 


V a 


C R P R’ V 


HA 


HC 


1 


1 


1 

K„-a 


h . • R . -a 

R 


A i 

C R ( V v - d 

+ 

R 1 

r n \r 




1 


l 

1 


h R’ a 

h * a 
+ c 


C R! , R -Vd 

c r p r‘ v 
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■m- 


87 


where : 

AV, AH 

QB 

HF 

HR 

HA 

HC 


v 

a 


d 

At 



h 


A 


h 


c 


h 


f 


HF 


K R-» + 
C R»R' v - d 


K, . a 
C A p A- Vd 


hj* a 

S 5 


effective heat 
el ements 
effective heat 
effective heat 
interface 
effective heat 
face 

effective heat 
air interface 
effective heat 
air interface 


conduction coef f iceints , neighboring 

generation coefficient, each element 
transfer coefficient, tire-rim 

transfer coefficient, tire-road inter- 
transfer coefficient, tire-outside 
transfer coefficient, tire-contained 


= volume of each element - see Table A 
= area across which heat is transferred - see Table A 
= distance between nodes - see Table A 
= time increment - selected 

= rate of heat generation each element - see earlier sections 

2 

= thermal conductivity of rubber = 5xl0 - ^ cal • cm/ cm -sec 

2 

= thermal conductivity of rim = .48 cal*cm/°C-cm -sec 
= specific heat of rubber = 0.. r cal/gm °C 

- specific heat of rim = .23 cal/gm °C 

3 

= density of rubber * 1 gm/cm 

3 

- density of rim = 2.7 gm/cm 

= film heat transfer coefficient, tire-road interface 

2 

- .25 cal/°C*cm -sec 

= film heat transfer coefficient, tire-outside air inter- 
face = ( 1 . 3x10^+4x10 ^V(cm/sec)) cal/ °C • cm^ • sec 

- film heat transfer coefficient, tire-contained air 
interface = 1x10 ^ cal/°C- cm^ • sec 

= film heat transfer coefficient, tire-rim interface 
= .5 cal/°C cm^-sec 
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Typical heat balance equations for the individual elements of 
the 22x5.5 8 PR tire are listed below where each equation has been 
solved for the change in temperature which takes place during a 
given increment of time. 

The inside elements all have an equation of the form: 

ATj i = QB( 1,1) + AHd.DCTj 2 -Tj j) + AV( 1 , 1 ) (T J+1 ^ ^ 1 ) 

+ AVO-l.lXTj^ j-Tj + HC(I)(TC-T I>;1 ) 

The middle elements all have an equation of the form: 

ATj g = QB( 1,2) + AH(I > 2 )(T i 3 -T i2 ) + AH(I, 1 )(Tj^-Tj 2 ) 

+ AV(I,2)(T i+ 12 -T i 2 ) + AV(I-1,2)(T i _ 12 -T i>2 ) 

The outside elements not in contact with the road surface or rim 
have an equation of the form: 

AT i>3 = QB( 1,3) + AH(I,2)(T i> 2 -T i>3 ) + AV(I,3 )(T i+ 13 -T i>3 ) 

+ av(i-i,3)(t i _ 1 > 3 -t i>3 ) + HA(I)(TA-T i>3 ) 

The outside elements in contact with the rim have the same form 
except HA is replaced by IIF and TA by TF. 

The outside elements in the tread region are assumed to be in 
contact with the air 90% of the time and in contact with the 
road surface 10% of the time. The heat balance equation for 
these elements is of the form: 

AT i 3 = QB ( I , 3 ) + AH(I,2)(Tj 2 -Tj 3 ) + AV( I . 3) (T y + 1 ^ g-Tj ^ g ) 

+ AV( 1-1,3) (Tj_ x 3 -Tj 3 ) + 0.9 HA(I)(TA-Tj 3 ) 


+ 0.1 HR( I ) (TR-Tj 3 ) 
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In this analysis the outside air temperature TA, the rim tempera- 
ture TF, and the road surface temperature TR are assumed to 
remain constant. However, the change in the contained air tempera 
ture TC is approximated as follows: 

h • A • At • (T. -TC) 

ATC = — 


where: 

C = .24 cal/gm ° c 

v 

p = .0013 gm/cm^ 

T. = weighted average of inside surface temperatures 

For the 22x5.5 8PR tire analyzed, the approximate values of A 

2 3 

and V are 30.4 cm and 90.35 cm 1 , respectively. Thus, for this 
tire: 

ATC = .1078(T -TC) • At 


